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A STOCHASTIC INVARIANTIZATION METHOD FOR ITO 
STOCHASTIC PERTURBATIONS OF DIFFERENTIAL 

EQUATIONS 

JACKY CRESSON, YASMINA KHELOUFI, AND KHADRA NACHI 


Abstract. In general, adding a stochastic perturbation to a differential 
equation possessing an invariant manifold destroys the invariance as far 
as the Ito formalism is used. In this article, we propose an invarianti- 
zation method for perturbations in the Ito case which can be used to 
restore invariance. We then apply our results to develop a stochastic ver¬ 
sion of the Landau-Lifshitz equation. We discuss in particular previous 
results obtained by Etore et al. in [6]. 


1. Introduction 

We consider a deterministic ordinary differential equation of the form 

(IT 

^ = f{x), hGN*. ( 1 . 1 ) 

A stochastic perturbation is taken into account by adding a ’’noise” term to 
the classical deterministic equation as follows: 

(iT 

- = /(x)+’’noise”. (1.2) 

and to replace the ’’noise” term by a stochastic one as 

dXt = f{Xt)dt + a{Xt, t)dWt, (1.3) 

where Wt is the standard Wiener process. This procedure is for example well 
discussed in ([10]). 

Of course, the main problem is in this case to find the form of the sto¬ 
chastic perturbation. We do not discuss this problem which is very compli¬ 
cated. We restrict our attention to the selection problem which is concerned 
with the characterization of the set of admissible stochastic models for a given 
phenomenon. By admissible we mean that the stochastic model satisfies some 
known constraints like positivity of some variables, conservation law, etc. This 
selection of a good candidate for a stochastic model of the phenomenon can be 
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done in many ways. However, in our particular setting, dealing with the sto¬ 
chastic extension of a known deterministic model, this selection is related to 
preserving some specihc constraints of the phenomenon. For example, part of 
the Hodgkin-Huxley model describes the dynamical behavior of concentrations 
which are typically variables which belongs to the interval [0,1]. This property is 
independent of the particular dynamics of the variables but is related to their 
intrinsic nature. The same is true for the total energy of a mechanical system. 
This quantity must be preserved independently of the dynamics. We formulate 
the stochastic persistence problem following our approach given in [5] in a 
different setting: 

Stochastic persistence problem : Assume that a classical ODE of the form 
(1.1) satisfies a set of properties V. Under which conditions does a stochastic 
perturbation of the form (1.3) satisfy also properties V ? 

The previous problem lead to characterize the set of a preserving the consid¬ 
ered properties V. classical properties are: Invariance of a given submanifold of 
MA, number of equilibrium points, stability properties of the equilibrium points, 
etc. 

The literature on invariance of manifolds for stochastic differential equations is 
huge and most of the time abstract in particular for what concerns the stochastic 
analogue of the Nagumo-Brezis theorem (see for example [2, 8]). This explain 
perhaps why these results are not so well known in the applied community be¬ 
cause the formulation of the conditions are not transparent for a given concrete 
system. In this article, we give a direct and simple derivation of a necessary and 
sufficient condition on the diffusion part in order that a submanifold globally 
defined as the preimage of a smooth function is preserved under a stochastic per¬ 
turbation. The result depends drastically on the stochastic differential framework 
that one uses. In the Stratonovich case, the condition on the diffusion is the same 
as the one on the drift part. However, in the Ito case, which covers most of the 
applications, the constraints for invariance are so strong that in many cases, one 
is unable to find a large class of admissible stochastic models. All these problems 
are discussed in Section 3. 

What to do if the framework for the model has to be the Ito one ? An idea 
is then to develop a systematic and algorithmic invariantization method in 
order to restore invariance in the Ito case. This is done in Section 4 following an 
idea initiated by Etore et al. [6]. In order to illustrate our method, we apply it 
to an Ito version of the Kubo stochastic Hamiltonian system and in Section 5.2 
to construct a stochastic version of the Landau-Lifshitz equation. We finally 
discuss the limitations and problems posed by the invariantization method. 


2. Reminder about Ito stochastic differential equations 

In this article, we consider a parameterized differential equation of the form 
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dXt = fit, Xt, b)dt, X G (DE) 

where b G is a set of parameters, / : M” x —> M” is a Lipschitz continuous 
function with respect to x for all b. We remind basic properties and definition of 
stochastic differential equations in the sense of ltd. We refer to the book [10] for 
more details. 

A stochastic differential equation is formally written (see [10],Chap.V) in dif¬ 
ferential form as 


dXt = fit, Xt)dt + ait, Xt)dWt, (IE) 

which corresponds to the stochastic integral equation 

fis,X,)ds+ f ais,X,)dWs, (2.1) 

Jo 

where the second integral is an Ito integral (see [10], Chap.Ill) and Wt is the 
classical Wiener process (see [10], Chap.II, p.7-8). 

An important tool to study solutions to stochastic differential equations is the 
multi-dimensional ltd formula (see [10], Chap. Ill, Theorem 4.6) which is stated 
as follows : 

We denote a vector of Ito processes by = iXt^i, Xt^ 2 , ■ ■ ■, Xt^n) and we 
put Wj = iWt,i,Wt, 2 , ■ ■ ■ ,Wt^i)to be a /-dimensional Wiener process (see [7], 
Definition 5.1, p. 72), dW^ = (dlTtq, dWt^ 2 , ■ ■ ■, dWt^i). We consider the multi¬ 
dimensional stochastic differential equation defined by (IE). Let E be a C^(M+ x 
M, M)-function and Xt a solution of the stochastic differential equation (IE). We 
have 

dF 1 

dFifXt) = -^dt + iX\F)dXt + -idXj)iXiF)dXt, ( 2 . 2 ) 

where VxE = dF/dX is the gradient of F w.r.t. X, V^E = VxV^E is the 
Hessian matrix of E w.r.t. X, 5 is the Kronecker symbol and the following rules 
of computation are used : dtdt = 0, dtdWt^i = 0, dWt^tdWtj = Sijdt. 

3. Stochastic invariance of snbmanifolds 

In this section, we derive an invariance criterion for a submanifold denoted 
by M of codimension 1 of MF which correspond to the zero set of a given function 
E : —)■ M of class i.e. 

M = {x G M”\E(x) = 0}, (3.1) 

under the flow of a stochastic differential equation in the Ito sense. This result is 
by itself not new and many general results are known in particular a stochastic 
Nagnno-Brezis Theorem as proved by Aubin-Da Prato in [2] or A. Milian in 
[8]. However, most of these results are difficult to read for a non-specialist in 
the field of stochastic calculus. The main interest of the following computations 
are precisely that our criterion can be easily derived using basic results in 
stochastic calculus. 
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3.1. Geometric definition of invariance. We consider an ordinary differen¬ 
tial equation of the form 


f xt = f(t,xt), 
\ x(0) = xo 


(ODE) 


where / : M"*" x —)■ M is a function of class and xq G is the initial 

condition. 


Definition 3.1. A given submanifold M C is said to be invariant under the 
flow of the differential equation (ODE) if for all xq G M, the maximal solution 
xt{xo) starting in xq when t = 0 satisfied xt{xo) G M for all t G M"*". 

We denote by TxM the tangent plane of M at x, we can write the invariance 
condition as follows 

f{t, x) G TxM, for all {t, x) G M"*" x M. (3-2) 

As M is of codimension 1, for all x G M we can define the normal vector N{x) 
to the tangent hyperplane TxM in x, such that 

TxM = {yGR^,y.Nix) = 0}, 

then the invariance condition can be written as 

N{x) ■ f{t,x) = 0, for all {t,x) G M"'' x M. (3-3) 

When M is of the form (3.1), the normal vector to M at x is equal to VE(x). 
Then the invariance condition reads as 

VF(x) • fit, x) = 0, for all (t, x) G M+ x M. (IF) 

In the stochastic case, the trajectories are continuous but nowhere differen¬ 
tiable. As a consequence the previous geometric condition can not be used. In 
the following we discuss two natural generalization of the notion of invariance in 
the stochastic setting. 

3.2. Strong stochastic invariance. Let us consider a stochastic differential 
equation of the form (IE). The stochastic character of the flow allows us to 
defined two natural notions of invariance. 

Definition 3.2 (Strong persistence). A submanifold M is invariant in the strong 
sense for the stochastic system (IE) if for every initial data xq G M almost surely, 
the corresponding solution x{t), satisfies 

F{F (x(t)) = 0,t G [to, +oo)} = I, 

i.e., the solution almost surely takes its values within the submanifold M. 

A direct computation gives the following criterion for stochastic invariance: 

Theorem 3.1 (Ito ’s strong invariance). Let M be a submanifold defined by a 
funetion F invariant under the deterministic flow associated to (DE), i.e., 

VF(x) • fit, x) = 0, for all X G M,t > 0 

The submanifold M is strongly invariant under the flow of the stochastic system 
(IE), if and only if, 

VF(x) • afl, x) = 0, for all x G M,t > 0 
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and 

Qx-dx ■ ^ 

i,j * 7 k=l 

Proof. The essential tool in this case is the ltd formula that will help us to for¬ 
mulate the invariance condition. Indeed, a process xt leaves the submanifold M 
invariant if and only if for all initial condition xq £ M a.s, the stochastic process 
associated to xt satisfies F{xt) = 0 for all t almost surely where it is dehned. 

The multidimensional ltd formula reads as 

^ Q'ijP 

d[F{xt)] = VF{xt)dxt + 2 X] Qy.,Q^ M t)dxi{t)dxj{f). 

i,j * ^ 

So we obtain 


d[F{xt] = VF{xt)f{t, xt)dt + VF{xt)a{t, xt)dWt+ 

+ 2 ^ dxidx ■ ^ 

i,j * ^ k=l 

The gradient of F being always normal to the tangent space of M, we have 
VF(xt) • f{t,xt) = 0 since the manifold M is assumed to be invariant in the 
deterministic case. It remains 

^ Q'ijP ^ 

d[F{xt] = VF{xt)a{t,xt)dWt + - ai^k{t,xt)aj^k{t,xt)dt. (3.5) 

i,j * ^ k=l 

The only contribution to the stochastic part is given by VF{xt)cr{t, xt) and is 
equal to zero if and only if the perturbation a satisfies the invariance condition 
(IF). Then the previous expression reduces to: 

^ q‘ 2 p ^ ^ 

d[F{xt\ = - ^ ^ Xt)aj^i{t, Xt)dt. (3.6) 

i,j * ^ 1=1 

that gives us the third condition. □ 

The previous Theorem indicates that unless a very specific form for a and F, 
there is no hope to recover invariance of a given manifold using a direct stochastic 
perturbation of a deterministic equation in the ltd case. 

As an example, we can specialize this result in the case of the sphere which 
will be important to study the invariance property of Landau-Lifshitz equation. 

Corollary 3.1. The sphere S is invariant under the flow of the stochastic system 
(IE) if and only if the stochastic perturbation is null on the sphere i.e., 

ai^i{t,x) = 0, i = 1, ...,n for all t£M.'^ and x £ S"'~^. 
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Proof. The proof follows from the fact that F(x) 
reduces to 


n 

xf so that condition 3.4 

i=l 


'^[(yi,i{t,xt)f = 0, V(t,x) e M+ X 5. (3.7) 

i=l 

This concludes the proof. □ 


As a consequence, trying to impose the invariance of S in the Ito case ’’anni¬ 
hilates” the perturbation that is intended to be produced by the diffusion term. 


4. The stochastic invariantization method 

In this section, we develop a procedure to restore invariance of manifold fol¬ 
lowing a procedure initiated by Etore and al [ 6 ] in a particular case. The basic 
idea is that in some cases, it is possible to project a flow which does not leave 
the manifold invariant on the manifold. 


4.1. First idea: projection procedure. Consider submanifolds of codimen¬ 
tion 1 of M”, that is defined by a homogeneous function of degree (7 G N;E : 
M" —)• M of class i.e., 

M = {xe W^/F{x) = 1} and F{\x) = \^F{x), for all x G A G M+. 

Let us assume that the coefficients of the system (IE) satisfy the invariance con¬ 
dition (IF), i.e., 

VF{x) ■ f{t, x) = VF{x) ■ a{t, x) = 0, for all x G M, t > 0, 
where a is a vector of MF and Wt is a scalaire Brownian motion, and assume that 


n 


E 


d^F 

dxidxj 


{xt)cFi{t, Xt)(7j{t, Xt)dt / 0. 


Then, by Theorem (3.1) we know that M is not invariant under the flow of (IE). 


A very simple way to construct an invariant stochastic process is to project on 
the manifold. In general, a projection on a manifold is difficult to compute. In 
our case, it reduces simply to consider the stochastic process 

yt = ^^, (4.1) 

F{xt)<i 

called the ’’projected” process associated to the stochastic process xt and the 
function F. 


Although simple, the previous method is in general not interesting. Indeed, 
the projected process satisfies in general a very complicated equation. 

Theorem 4.1. Let us assume that xt is the solution of ltd equationiJE) and M is 
defined by F of class as above. The projected process yt satisfies the equation 
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dyk = 


F{x) 


1 dF 
2q dxk 


. 1+2 


(x).F(x) 2 (Tk{t,x) 




1 -y ^ — 1 

+ 9 X] —Xkcri{t,x)aj{t,x) 




_i+2 d'^F 


dxidxj 


1 + g l+^dF dF 

-F(x) 2 —(x) —(x) 

q oxi oxj 


dtF 


+ 


F{x) 2 cjfc(t,x) 


dWt, for each k = 1, 


We leave the proof to the reader as it follows easily from the ltd formula. 

Even for a simple manifold as a sphere, the corresponding equation does not 
simplify. 

Corollary 4.1. Let us assume that xt is the solution of ltd equation (IE). The 
projected process yt = —on the sphere satisfies the equation 

F{xt)^ 


dyk = 


F{x) 2/fc(t,x) - ^XfcF(x) 2crfc(t,x)+ 


+ E ai{t,x)aj{t,x)F{x) ^XiXj 




dt 


(4.2) 


+ 


F(x) 2 ak{t,x) 


dWt] for each k = 1, 


We leave the proof to the reader. 


The previous expression has many problems: 

• The resulting stochastic differential is far from being simple and can not 
in general be written only with respect to yt- 

• The form of the deterministic part can not be seen as a perturbation of 
/. This induces difficulties for the interpretation of the new equation. 

In the following, we follow a different strategy initiated by P.Etore et al. in 

[ 6 ]. 

4.2. The invariantization method. We first introduce the notion of invaria- 
tized process. 

Definition 4.1 (Invariantized process). Let xt be a diffusion process defined by 

dxt = f{t, xt)dt + a{t, xt)dWt. (4.3) 

The invariantized process associated to (4.3) and the submanifold M defined by 
F is defined by 

f dyt = f{t, xt)dt + a{t, xt)dWt 

Xt = —(4 4 ) 

(Fiyt)F ^ ’ 

yQ = y e M 
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This terminology is justified by the fact that we have 

F{xt) = 1, for all f > 0. 

The method associating to a given process and a submanifold M its invari- 
antized process is called the invariantization method. 

The main property of the invariantization method is that the stochastic differ¬ 
ential equation satisfied by xt is simple in the contrary to the projection method. 

Theorem 4.2 (Invariantization). Assume that F is homogeneous of degree q and 
that F{yt) is a non random process, i.e. that there exists a function h{t) such 
that 

dF{yt) = h{t)dt. (4.5) 

Let us denote by H (t) the function defined by 

t 

(4.6) 

0 

Let us assume that L[{t) >0 for t > 0, then the invariantized stochastic process 
associated to yt and F satisfies the stochastic differential equation 


L 

H{t) = 1 + J h{s) ds. 


dxt = 


ifm 

qH{t) 


Xt + 


1 


iH{t)y 


-f(t,xt) 


dt -\- 


1 


iHit)y 


—a{t,xt)dWt 


(4.7) 


where H denotes the derivative of H with respect to t. 


Proof. This is a simple computation. As F{yt) = H{t), the stochastic process 
Xt = yt/H{ty/^ satisfies 


dxt = - - 


1 H{t) 


-yt + 


-dyt- 


qH{ty/<i+^^ H{ty/<1 

Replacing yt by xtHity/^^ and dyt by its expression, we obtain 

+ -^—Tfit,xt) 


(4.8) 


dxt = 


This concludes the proof. 


H- ^-^a{t,xt)dWt. 


dt 


□ 


5. Applications: an Ito Kubo oscillator and a stochastic 
Landau-Lifshitz equation 

5.1. Example: Ito version of the Kubo oscillator model. Let us consider 
the Kubo oscillator (see for example [9]) in the Ito case, which can be written as 

dXt = JaXtdt + JMWt, (5.1) 

where X = G M^, a, ct G M, Wt is a 1-dimensional Brownian motion and 

*=(" 
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The Stratonovich version of the Kubo oscillator has any circle Xf + 

Vr G M"*", invariant under the flow. However, the circles are not invariant under 
the flow of the ltd version of the Kubo oscillator. Indeed, we have using the ltd 
formula with T’(Xi, X 2 ) = Xf + X|, that 

dF((Xi, X 2 ) = a‘^{Xl + Xl)dt. (5.2) 

Assuming that T’(Ai, X 2 ) = r^, r 7 ^ 0 is invariant under the flow gives 

dF{{Xi,X2) = a^r^dt. (5.3) 

As a consequence, the condition dF = 0 is satisfied if and only if u = 0. This 
means that in the ltd case, invariance can not be preserved while the flow is 
stochastic (i.e., cr / 0 ). 


If we apply the last transformation such that Xt = ppr we find that Xt is the 


solution of the stochastic system 


dXt — Ja,a,t Xfdt + 


\/ + 1 


\Yt\ 

JaXtdWt, 


(5.4) 


where 


— 1 g __g2_ 

\ Vcr^t+1 2(cr^t+l) ^ 


which preserve the invariance of S under the flow of the deterministic equation 

dXt = JaXtdt, for all a G M. (5-5) 


5.2. An ltd stochastic Landau-Lifshitz equations. In this Section, we de¬ 
rive a stochastic Landau-Lifshitz equation in the ltd setting. We first remind the 
construction of the classical Landau-Lifshitz equation and then its main proper¬ 
ties. We then review classical stochastic approach used by different authors and 
the difficulties associated with these models. 


5.2.1. The Landau-Lifshitz equation. The Landau-Lifshitz equation is a general¬ 
ization of the classical Larmor equation. The Larmor equation is conservative. 
However, dissipative processes take place within dynamic magnetization 
processes. The microscopic nature of this dissipation is still not clear and is 
currently the focus of considerable research [1, 3]. The approach followed by 
Landau and Lifshitz consists of introducing dissipation in a phenomenological 
way. They introduce an additional torque term that pushes magnetization in the 
direction of the effective field. The Landau-Lifshitz equation becomes 

= —ji X b — afj, X (ij, X b), (LLg) 

where ^ G is the single magnetic moment, x is the cross product in b is 
the effective field and a > 0 is the damping effects. 


As for the Larmor equation, this equation possess many particular properties 
which can be used to derive a stochastic analogue. We review some of them in 
the next Section. 
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5.2.2. Properties of the Landau-Lifshitz equation. In this Section, we give a self- 
contained presentation of some classical features of the LL equation. Readers 
which are familiar with the LL equation can switch this Section. 

Invariance. The following result is fundamental is all the stochastic generalization 
of the LL equation. 


Lemma 5.1. Let /i(0) G 5^, then the solution pit satisfies for all t G 
i.e. the sphere is invariant under the flow of the LL equation. 


ht 


1 = 1 , 


(5.6) 


We give the proof for the convenience of the reader. 

Proof. Let nt be a solution of the LL equation. We have 

= Ht. [-Ht X 6 - ant X (nt X 6)]. 

By definition of the cross product, the vectors fit x b and a fit x {fit x b) are 

orthogonal to fit so that 

^ [fit.fit] = 0. (5.7) 

As a consequence, using the fact that fio £ S'^, we deduce that 

WhtW = ll/^oll = 1, (5.8) 

which concludes the proof. □ 

As a consequence, a solution starting on the sphere will remains always on it. 
The sphere being a two dimensional compact manifold, we can use classical result 
to deduce the asymptotic behavior of the solutions. But first, let us compute the 
equilibrium points. 

Equilibrium points. The equilibrium points of the LL equation are easily ob¬ 
tained. 


Lemma 5.2. The LL equation possesses as equilibrium points 6/ 

-b/ II b ||. 


and 


We give the proof for the convenience of the reader. 

Proof. An equilibrium point /x G satisfies 

—fi X b — afi X {fi X b) = 0, (5-9) 

which gives 

—fi X b = afi X {fi X b). (5.10) 

The vector fi must be orthogonal to fi x b and at the same time equal to —fi x b 
up to a factor a > 0. As /xq G 5^, we have fi 0 and the only solution is 

fixb = 0. (5.11) 

We then obtain fi = Xb, with A G M. By Lemma 5.1, we must have fi £ S'^ so 
that A = ±1/ II 6 II. This concludes the proof. □ 

We see that the equilibrium point of the LL equation coincide with those of 
the Larmor equation. 

The stability of the previous equilibrium point can be easily studied using the 
Lyapunov theory. 
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5.2.3. Toward a stochastic Landau-Lifshitz equation. In this Section, we discuss 
the usual way of deriving a stochastic analogue of the Landau-Lifshitz equation 
by considering an external perturbation of the effective magnetic field. We focus 
on the Stratonovich and the Ito interpretation and we explain the strategy used 
in Etore et al. [ 6 ] to bypass the obstruction that the Ito version does not preserve 
the sphere using the invariantization method. 

Classical approach to the stochastic Landau-Lifshitz equation. The main ap¬ 
proach to deal with the stochastic behavior of the effective magnetic field is to as¬ 
sume that the effective field b is subject to a stochastic perturbation b+ €” noise”. 
Due to the linearity with respect to the parameter b, we obtain an equation of 
the form 

d/Jt = [—ht X b — afit X {lu-t X b)]dt-\-£[—/at x ”noise” —aixt x htX ”noise”]. (5.12) 

Interpreting the previous equation in the Ito formalism of stochastic differential 
equation leads to the following stochastic model; 

dqit = [—Tt X 6 — aqit x (/Uf x b)]dt -|- e[—/r* x dWt — ant x x dWt], (ELL) 

where the term cj(t, x) = —x x • — ax x (x x •) can be written as 

( a(x| -|- x|) X 3 — axiX 2 —X 2 — 0 x 3 x 1 

—X 3 — 0 x 1 X 2 o(x| -|- x\) xi — 0 x 3 X 2 

X 2 — 0 x 1 X 3 —xi — 0 x 3 X 2 o(x| -|- xf) 

Most authors use the Stratonovich formalism in order to give a sense to the 
previous equation. The main reason is that in this case, the invariance of 5^ is 
ensured. However, as pointed out by Etore et al. in [6], the Stratonovich ver¬ 
sion of the Landau-Lifshitz equation leads to several difficulties, such as the fact 
that the stability of the equilibrium points of the deterministic LL equation is lost. 

The previous point has motivated the work [6] in which the authors discuss the 
Ito case. However, the Ito approach lead to other difficulties. Details are given 
in the next Section. 

Stochastic Ito perturbation of the Landau-Lifshitz equation. The Ito version of 
the stochastic Landau-Lifchitz equation possesses many drawback and the main 
one follows directly from the invariance criterion that we derive in Corollary 3.1. 
Indeed, we have: 

Lemma 5.3. The sphere 5^ is not invariant under the flow of the ltd version of 
equation (ELL). 

Proof. By Corollary 3.1, the diffusion term must be zero on the sphere S'^ and 
all t £ M. The condition on a on the diagonal terms implies that a = 0, i.e. 
that we can not have a dissipative term and we recover the Larmor equation in 
contradiction with our assumption that a 0. This concludes the proof. □ 

The previous result excludes the use of the Ito formalism for a direct stochastic 
generalization of the Landau-Lifshitz equation. However, we can use the invari¬ 
antization method exposed in Section 4 in order to obtain an Ito model, related 
to the previous one, but which satisfies the invariance of the sphere S'^. 


(5.13) 
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Invariantization of the Landau-Lifshitz stochastic Ito model. Let us consider the 
stochastic LL equation (ELL). The sphere 5^ is defined by the homogeneous 
function F{xi,X 2 ,x^) = degree 2. Let us consider the invariantized 

process defined by 

J dyt = [-yt xb- ayt x (fj-t x b)]dt + x dWt - afit x //* x dWt], 

I = ik- 

(5.14) 

A simple computation gives (see also ([6], Proposition 1)): 

Lemma 5.4. The process F{yt) is random. Precisely, we have dF{yt) = 2e^(a^ + 
l)dt. 


As a consequence, Theorem 4.2 applies and we have: 


Lemma 5.5. The invariantized stochastic differential equation associated to 
= Ylkik ^be ltd stochastic differential equation (ELL) is given by 


dxt = 


1 2e^(a^ + 1) 1 

-i— Xi A - - 

2 2e2(a2 + i)t + i ^2e2(a2 + i)t + i 

H — / r,, } , tj(f, xt)dWt, 

V2e2(a2 + i)t + i 




dt+ 


(5.15) 


Proof. This is a simple computation. □ 

We recover the Etore et al. version of the Stochastic Landau-Lifshitz 
equation proposed in [6]. 


5.2.4. About equilibrium points. A natural question about the previous invari¬ 
antized model is to up to which extent it answers to the reasonable constraints 
one waits for a stochastic version of the Landau-Lifshitz equation. Eor example, if 
one is interested in preserving equilibrium points of the initial system, the model 
is not satisfying. Indeed, we have: 


Lemma 5.6. The points fi = ±6/||6|| are not equilibrium points of equation 
(5.15). 

Proof. Eor fi = ±6/ || b ||, we have for all u G that 

a{t, ±6/ II b II).u = “11““))^ XV — ±a||-^ x (±||-^ x v). (5.16) 

The second term is always zero but the first one is only zero when v is co-linear 
with b. However, as v takes arbitrary values, we can not ensure this equality. As 
a consequence, the initial equilibrium points are destroyed under the stochastic 
perturbation. □ 


It must be noted that the previous problem can be easily solved by modifying 
a little bit the modeling of the stochastic behavior of the effective magnetic held. 
Indeed, let us consider instead of dWt the following vector 

bdWt, (5.17) 

where Wt is a one dimensional Brownian motion. This assumptions is equivalent 
to say that we consider only stochastic behavior in the direction of the initial 
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field b. This is of course very particular, but in this case the new model preserves 
the equilibrium points of the initial system: 


Lemma 5.7. Let us consider the modified Etore et al. stochastic Landau-Lifshitz 
equation defined by 


+ 


+ l)t + 1 


dxt = 




1 2e2(a2 + 1) 

~2 2e2(a2 + l)t + 1 
dt -\- 


Xt + 


^2e^{a^ + l)t + 1 


a{t,xt) ■ [bdWt] , (5.18) 


where a is defined by equation (5.13) and Wt is a one dimensional Brownian 
motion. This equation possesses as equilibrium points ±6/||6||. 


Proof. This follows easily from the previous proof only saying that v is always 
co-linear to b. □ 


6. The invariantization method as a stochastic perturbation 

Although the invariantization method leads to a simpler equation than the pro¬ 
jection procedure, it is not very easy to understand the procedure as a stochastic 
perturbation of the deterministic model. In this Section, assuming that the dif¬ 
fusion is governed by a small parameter 0 < e << 1, we write the invariantized 
equation as a perturbation. 

6.1. Small perturbation and invariantization. In the following, we use the 

notations of Section 4. For <7 = 0, the invariantized process reduces to the 
deterministic equation. Let us assume that a is of the form 

a{x) = ecro(x), 0 < e << 1, (6.1) 

where do and / satisfy the invariance conditions. 

Using the Ito formula, we obtain 

dF{yt) = VF{yt).dyt + ^d'^Fdy'^dyt.dyt. (6.2) 

As dyt.dyt = e‘^aQ{xt)dt, and F and do satisfy the invariance relation, we finally 
obtain 

1 d^F 

dF{yt) = {xt)dt. (6.3) 

Denoting by 7 ( 1 ) the function 

t 

1 r n2 p 

lit )=2 y -^^ii^t)dt, ( 6 . 4 ) 

0 

we then obtain using Theorem 4.2 a function h^ of the form 

hfit) = e^-f{t), (6.5) 

and a function Fl^ of the form 

(6.6) 


Hfit) = l + e^d{t). 
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As a consequence, we can develop the drift part with respect to e and we obtain 
for the invariantized process an equation of the form 

dxt = (j" + * + . ..^dt+ ^ecTo + + ... ^ dWt- (6.7) 

We do not search for explicit expression of the perturbation terms * and *. 

6.2. Limitation of the method. From a modeling point of view, we believe 
that a stochastic model of a deterministic equation must satisfy the following 
constraints: 

• First, the drift part must be relied early to the deterministic equation and 
moreover must be understandable as a perturbation of it, i.e. of the form 

f{t,x{t)) + P{t,x{t)) (6.8) 

where P{t,x{t)) is the perturbation term. 

• Second, the new equation must be easy to interpret and must keep a sense 
with respect to the field of applications. 

What can be said about the invariantization method ? 

The first modeling constraint is then satisfied by our invariantization method. 
However, the perturbation term obtained in equation (6.7) is very complicated 
and the role of each term in the dynamics is not easily recovered. 


7. Conclusion and perspectives 

For ltd stochastic perturbation of ordinary differential equations, we have de¬ 
rived a general method allowing to preserve invariance of a given codimension 
one submanifold under the stochastic flow. This method has however some lim¬ 
itations and lead to difficulties in the interpretation of the resulting model from 
a perturbative point of view. A natural question is then to find other stochas- 
tization procedure which still use the ltd formalism for stochastic differential 
equations but for which invariance can be ensured under reasonable constraints. 
We refer to [4] where this problem is discussed in general in the framework of 
random ordinary differential equations and used to construct a new model of a 
stochastic Landau-Lifshitz equation. 
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